^ 000659446 >ns. vol. -ts. no. y may im 



p.d. Aj|bi83L. , , 




A Novel Block Truncation Coding of Color Images 
Using a Quaternion-Moment-Preserving Princip^ 



Soo-Chang Pei, Semor Member. IEEE, and Ching-Min Cheng, Member, IEEE 



Abstract — Block truncation coding (BTC) is an efficient tool 
for image compression. To compress color-pixel blocks, a novel 
color BTC algorithm, called quaternion-moment block truncation 
coding (QMBTC), is presented in this paper. Analytical formulas 
for QMBTC, whose computation time is on the order of pixel 
block size, are derived by using quaternion arithmetic and the 
moment-preserving principle. The proposed color BTC algorithm 
can adaptively truncate a pixel block into one or two output 
classes according to the distribution of color values inside the 
blocks. The experiments] results show that the compression ratio 
is increased as compared with existing color BTC algorithms, 
and the picture quality of reconstructed images is satisfactory. 
In addition, a post-BTC data compression scheme is proposed 
to further compress the sub image constructed by reproduction 
colors of truncated pixel blocks. Using a lookup table to display 
decoded data, this postprocessing scheme can output images 
acceptable to human eyes. 

Index Terms — Block truncation coding, color image compres- 
sion, quaternion-moment-preserving. 



I. Introduction 

BLOCK truncation coding (BTC) was first proposed by 
Delp and Mitchell [1], [2] to compress monochrome 
images. Unlike other image compression methods such as 
transform coding and vector quantization [3], BTC requires 
less computation effort and has good capability of combating 
channel errors. Lema and Mitchell also proposed a mono- 
chrome BTC algorithm [4] by preserving the mean-absolute 
moment or minimizing the mean-absolute error in each pixel 
block. The outputs of these BTC quantizers for each pixel 
block include two numbers which are the associated centroids 
of two partitioned classes and a bit map which specifies the 
states of each output pixel. The bit map (1 bit/pixel) occupies 
more than 61% of the output code. Although some efforts have 
been made to compress the bit map [5], [6], the subjective 
results are not quite satisfactory. 

The BTC method described above can be extended to color 
images by applying the BTC technique to each of the color 
planes [4]. Based on empirical results [4], good quality images 
at 6 bits/pixel have been reported. Since BTC is a two-level 
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quantizer that adapts to local properties of the image, the tHree 
resulting bit maps will be quite similar or almost identical. 
This motivates the use of one bit map to quantize all three of 
the color planes. Thus, only one out of three bit maps needs 
to be stored. Several single bit-map color BTC algorithms 
[7H9] have been proposed. The color cell compression (CCC) 
algorithm, which is presented by Campbell ex at. [7], uses the 
sample mean value of the luminance values of the pixel block 
to truncate the color pixels inside the block into two classes. 
Then the centroids of each class are chosen as the reproduction 
colors of the block. Wu and Coil's algorithm [8] truncates the 
average of the R, G, B planes for the corresponding pixels by 
choosing the sample mean value as the threshold. Recently, 
Kurita and Otsu proposed a color BTC algorithm [9], which 
is based on truncating the principal component of R, G. and 
B values inside the pixel block. They reported that the mean 
vector is preserved, and the covariance matrix of the decoded 
block is the same as the between class covariance matrix. 
When no post-BTC data compression scheme is employed, 
these techniques [7]-[9] can obtain the resultant bit rate of 4 
bits/pixel for the 4 x 4 block size. 

In contrast to the above-mentioned color BTC algorithms, 
which process on only one transformed component of input 
color data, we propose in this paper a novel color BTC 
algorithm, called quaternion-moment block truncation coding 
(QMBTC). QMBTC generalizes conventional monochrome 
BTC [ 1 ] to color BTC by expressing the input color space as a 
quaternion-valued space. Through the definition of quaternion 
moments of input color data, we extend the moment-preserving 
principle introduced by Delp and Mitchell [1] from one- 
dimensional (1-D) monochrome data to three-dimensional 
(3-D) color data. An analytic solution for QMBTC is also 
obtained by the usage of quaternion arithmetic. One feature of 
QMBTC is that it can determine the number of output classes 
according to the distribution of color values inside the pixel 
block. A pixel block with similar color values produces only 
one output class, and the associated bit map can be replaced 
by a I bit reference indicating that one color clustering 
happened. Thus, QMBTC can achieve a better compression 
ratio than existing color BTC algorithms. As well, with no 
severe performance degradation, a post-BTC data compression 
method is introduced to further compress color images. 

This paper is organized as follows. Some preliminaries of 
quaternion arithmetic and the definition of quaternion moments 
are first introduced in Section II. Section III presents the 
QMBTC and its performance evaluation. In Section IV, the 
application of the proposed single bit map color BTC for image 
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compression is addressed, including the implementation of a 
post-BTC data compression method. In Section V, we compare 
the performance of existing color BTC algorithms with the 
proposed algorithm. Finally, some conclusions are made in 
Section VI. 

II. Quaternion Moments 

Human perception of gray level images is quite satisfactory, 
but color images do seem to be perceptually richer. And, 
different from gray level, color data representations are usually 
ternary. For example, color image data from a frame grabber 
are usually in the RGB (red, green, blue) primary color space. 
There are other colorimetric representations which allow the 
entire visual gamut being described. Some of the examples 
are: C1E1930-XYZ, CIE1976-L*u*6\ HSV, and CCIR601- 
YCrCb color coordinate spaces [3]. To effectively compute 
multidimensional color data, an efficient expression of color 
data is necessary. Machuca and Phillips [10] have proposed 
vector fields as the theoretical model of color data. Through 
the usage of differential geometry and vector analysis, they 
obtained the algorithm for color edge detection. 

In contrast to using vector fields [10], we select an algebraic 
approach in this paper to generalize the moment-preserving 
principle of monochrome BTC [1] to multidimensional color 
image. Although real and complex number systems are used 
to provide arithmetic operations of one-dimensional and two- 
dimensional data, there are no usual laws of arithmetic to 
handle algebraic operations of ternary numbers. We choose 
the quaternion number, which was discovered by Hamilton 
[II], as the expression of color data. In this expression, one 
dimension of the quaternion number is redundant and set 
to zero. The algebra of the quaternion number, quaternion 
arithmetic, is the generalization of complex numbers. Using 
quaternion arithmetic, analytical formulas for monochrome 
BTC [1] can still be maintained for a color image. Routine 
BTC [1] is a special case of the proposed moment-preserving 
BTC in the quaternion space. In this section, we will address 
some preliminaries of quaternion arithmetic and the definition 
of quaternion moments used by the following sections. 

Considering a 4-D real-valued data set H = 
ffi(»).92(n) 9 to(»))}»«i- a quadruple data point 
(9o("Mi(rc),<72(n),g3(«)) can be expressed as a quaternion 
number q(n) 

q(n) = q 0 (n) + 7i(n) • < + q2(n) • j + 7a(n) • k (I) 

where i.j, and it denote the operation units of quaternion 
number. Any vector v € R? can be expressed as a quaternion 
with g 0 set to be zero. For example, a color value (R.G,B) 
can be shown as a quaternion with q% = R*q 2 = G, 73 = 
B,qo = 0. Any vector v € R 2 can be expressed like a complex 
number. A quaternion can also be denoted as q(n) = (0.6) 
where a = {qi(n),q 2 {7i).qz{n)) and * = Qo( n )- 

Operations on quaternion number have the following prop- 
erties. 

1) The addition and subtraction rules of quaternions are the 
same as for complex numbers. 



IFJ-JS TRANSACTIONS ON COMMUNICATIONS. VOL. 4.V NO. 3. MAY .|««f 

i. * '"• * 

2) Using the cross product of vector space ( x ), one can 
define multiplication of two quaternions, ij and 7' as 

q<<j'=(aJ>)-(a'.h') 

= <a x a' + /j • a' + b' ■ a. h - // - a • a'). (2) 

3) The conjugate //* of q is defined as 

t) m - -(a. h) = f/ n - (q\ ■ i + q 2 • j + 7.1 " *•) < 3 > 

and the norm of the quaternion is denoted as ||#/|| 2 = 7 • 7*- 

4) The reciprocal of q is 

With the help of the (7)" 1 , the division of quaternions is 
denoted as 

i=q'(qr l . (5) 
Q 

Based on the above definition of the quaternion, we will desig- 
nate the first-, second-, and third-order quaternion moments as 
follows in order to explicitly express the statistical parameters 
of a 4-D data point: 

mi = E[q] 
m 2 = E[q - q*] 

r?*3 = E[q -q* q] (6) 

where E[-] represents the expectation. The definitions of m\ 
and rh2 are me extension of complex moments. The definition 
of the third-order quaternion moment rhs is adopted from 
high-order statistics [12]. Equation(6) can be further expressed 
as 

mi = E[q 0 ] + E[q{\ • * + E[q 2 ] • j + Efo] ■ k 

m 2 = E[q$ + q\ + 75 + q\] 

m 3 = E[ql + q\qo + 7?7o + <?a<?o] 

+ E[q\ + qlq x + q^qi + q%qi] • i 
+ E[q\ + qlq 2 + q\q 2 + ql<n\ ' j 
+ E[q$ + 7573 + <?i<?3 + q 2 qz\ • k. (7) 

To calculate the quaternion moments of the data set H, the 
expectation !?(•] of (7) is replaced by the sample mean value. 
As we observe, moment m L represents the average magnitude 
of the input data (70* <7i ■ ?2« 9a)> The real- valued m 2 expresses 
the expected value of the power of (goitfi^^tfa)- The third 
moment m a consists of the sum of joint third-order moments 
among 7o« <7i • '/2< an d 73 which reflect the measure of skewness 
of the input data. Although statistical distribution of a color 
image can be described by knowing each order of moment, 
it is known that the higher the order of moment is, the less 
important this moment is. Thus, preserving moments up to 
third order would not severely lose statistical distribution of 
the image. To derive analytical formulas for QMBTC in the 
next section, we will use the moment-preserving principle on 
m L ,m 2 i and mj in order to maintain statistical characteristics 
of input image. As empirical results will show, the output 
picture quality is satisfactory for human eyes. 
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J III. Quaternion-Moment Block Truncation Coding 

The problem of applying QMBTC in a pixel block H is 
to sel ct a hyperplane as a threshold such that, if all below- 
threshold pixels in H are replaced by the quaternion value 
i 0 and those above-threshold pixels in H are set to be the 
quaternion value £i, then the first three quaternion moments of 
image H are preserved in the resultant two-level pixel block G. 

Let po and p x denote the probabilities of the below-threshold 
and the above-threshold pixels in H after the QMBTC, respec- 
tively, and 



Po + Pi = 1. 



(8) 



Then the first three quaternion moments of the two-level pixel 
block G are given by 

l 



i=0 

1 



TO3 



1 

> , Pi ' *\;*\ ■ zi- 



(9) 



1=0 



If we let the first three quaternion moments of the two- level 
pixel block G be equal to those of H, 771 1, = 771 J,. U = 1,2,3. 
we get the following moment-preserving equations: 



PO • ZQ +Pl • ii =771! 
PO • ^0 • ^0 + Pi ' *1 * *l = ™2 
*0 ' ZQ * *0 + Pi * *1 • ij * Zi = 77i 3 . 



(10) 



Thus, using the quaternion moments, the Delp and Mitchell 
moment-preserving principle [1] can still be maintained for 
color image data, and the solution will be a special case of 
the proposed QMBTC. 

The moment-preserving equations (10) can be solved in- 
directly by adopting the following polynomial of quaternion 
z. 

C(z)=.z z m + z* -Ci + Co 

= (i-£o)(£-ii) (11) 

where i; for / = 0, 1, the quantized levels of QMBTC, are the 
roots of C(£), i.e., C(zi) = 0. If we multiply C(z t ) with (8) 
and the equation defining rhj, respectively, we get 

\po • £0 ' C(i 0 ) +Pi • zi • = 0. 

Then, using (10), Co and ci of (11) can be obtained as 

_ (77I3 — TTli ' ^2) 

(rh\ • 771 j — 7712) 
co=-(m;.ci+7h 2 ). (13) 



(12) 



ci = 



With the help of Co and ci, the roots of C(£/) = 0 can be 
solved. 

Considering the case of thresholding color pixels with q 0 
set to zero, the unknown zi can be denoted as 



Zi = Zll ■ i + Zf2 * i + 2/3 • *. 



(14) 



If we express ci and r 0 of (13) as 



f'o = <*w) + t'tn ■ i + 012 ' i + ' : o;i * k 



(15) 



we have the following results after putting (14) and (15) into 
(It) and letting (II) equal zero. 

^12 * z\z + cot 
~n = 

r Z2 = 

a * zj^ + b • zrj + d = 0 



(16) 



where 



rt=rj A +c72+Cuj 



6 = 2 • cox • C12 — 2 • C02 * Cn + c 2 x • r l3 

c/ ss r HU • C12 * Cu — C02 ' Cu • CJ3 
+ C 01 + C00 -cf 3 . 

From (16), the two solutions z 0 3 and z i3 can be obtained 
directly from the quadratic formula 



-1 



z l3 = (6 + \/6 2 -4-a-ei). 

J ■ a 



(17) 



The corresponding zn and 212 can thus be acquired from (16). 
After £ 0 and z x are obtained, the h yperplane which is 
perpendicular to the line segment £ 0 ii such that the quaternion 
space is split into two halves, can be defined as the decision 
boundary to separate the input data set into two classes. The 
half space containing £ 0 has po portion of pixels, and the other 
half space containing £ t has p x portion of pixels. However, if 
we want to implement the QMBTC in a much quicker way, 
/' can be chosen alternatively as a hyp erplane perpendicular 
to and bisecting the line segment £ 0 £i» If the color value of 
each pixel is denoted as (h, J 2 ,, JaJ/the linear. equation, which 
formulates the decision boundary is thus described as 



I\ = s — t\ • I2 — *2 * J3 



(18) 



with 



8 = 



1^2 



± 4> + gff a - z n - 



12 



ti = 



(^02 ~ ^12) 
(201 - :u) 

* ( r 03 ~ ^13) 

h = 7 r. 

(-01 - z n ) 

This segmentation of two classes is equivalent to the nearest 
neighbor clustering of two centers i 0 and £ x . The classes thus 
obtained can be used to produce a single bit map for the pixel 
block. Besides, the computation complexity of the QMBTC is 
mainly dominated by the calculation of quaternion moments 
7h„. for v = 1,2,3. whose computation time is of order /V, 
the size of pixel block. 
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A. The Algorithm 

In summary, the proposed color BTC algorithm can be 
described as follows. 

1) Divide the input color image into small nonoverlapping 
blocks, 4x4 size in this paper. 

2) Express the color value of each pixel (A. /a, h) by the 
quaternion number 

q = <7o + 9i ' » + 92 i + <7:j ' * 

with tfi = /i , <7 2 = /a. 9a = h,Qo = 0; i,j, and ^ denote the 
operation units of the quaternion number. 

3) Obtain two quantized levels i 0 and fx for each pixel 
block by solving the quaternion-moment-preserving equations 
(10). 

4) Choose the hyperplane /' defined by ( 1 8) as the decision 
boundary to segment the pixel block into two halves. 

5) Construct a class-indicating bit map such that each pixel 
location is coded as a **one" or a "zero," depending on whether 
or not the pixel resides in the half containing Si. 

Instead of using zq and ii, we select in this paper the 
centroid of each output class as reproduction colors of the trun- 
cated block in order to reduce the minimum mean-square error. 
To illustrate coding a pixel block, we select a 4 x 4 pixel block 
from an image and arrange the color values (/i.Ja,^) = 
(R,GiB) in the block as a quaternion-valued matrix X. In 
X. each element (q 0 , q\ , $2* Qz) is set to (0. R< G,B) of the 
corresponding pixel location, as shown on the next page. 

i r = (0, 229.6, 135.5, 109.2) 
i 0 =(0,206.1,87.6,9.4) 

and the decision boundary /' is h = 697 - 2.03 • J 2 - 4.24 - J 3 . 
By deciding whether or not the value of I\ - 697 + 2.03 • 
/ 2 + 4.24 • / 3 is greater than zero for each pixel, we have the 
corresponding bit map expressed as 

\i i i i/ 

The root mean square error (RMSE) of QMBTC between the 
original and quantized block can be calculated to have the 
value 7.71. This example shows that only one output class 
is generated. When employing a A:-means classifier [13] to 
truncate the same pixel block, we arbitrarily choose two pixel 
values as the initial values. After three iterations, we obtain 
two output classes which can be illustrated by the bit map 

(1 1 0 0\ 
0 1 0 0 
0 0 1 1 I ' 
0 0 0 1/ 

If each pixel inside the corresponding class is represented by 
its centroid, the RMSE of the A:-means classifier has value 
6.03. Another example X\ shown on the next page, will 
demonstrate that QMBTC produces two output classes: 
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ii =(0., 146.0,85.2, 102.2) 
£ n =(0., 14.1. 13.3.29.3) 

and the decision boundary /' is /i = 143.2 — 0.545 •/;> — 0.552 • 
/a with the corresponding bit map expressed as 



\0 1 1 1/ 

The RMSE of QMBTC has value 14.02. Using the Ar-means 
classifier [13] in this example, we also arbitrarily choose two 
pixel values as the initial values. After five iterations, we 
obtain the same two output classes as those of QMBTC. 

From these examples, we observe that if color values of 
the block are similar, the QMBTC will return one output 
class. Otherwise, if the color values of the block are not very 
close, two color clusterings will be obtained. Concerning the 
A:-means classifier, it always produces two output classes and 
the complexity of the A:-means classifier is greater than that 
of QMBTC. For instance, the A:-means classifier should worry 
about the number of iterations needed and the problem of being 
convergent or not. For the data size and the number of iter- 
ations being equal to N and /, respectively, the computation 
load of the ^-means classifier is of order J • iV*. In addition, 
unless the initial guess is well chosen, the results of the fc- 
means classifier are easily trapped inside the local minimum, 
which degrades the performance of the fc-means classifier. 

If the A-means classifier wants to adaptively truncate a pixel 
block into one or two output classes, one intuitive measure 
is the RMSE between the original and quantized block. As 
we see, RMSE's of the fc-means classifier for the above two 
examples are not different significantly. Moreover, considering 
a one-class block, whose pixel values are near, and a two- 
class block, whose pixel values are equal for each class, we 
can expect that the RMSE's of the fc-means classifier for these 
blocks are very close. Thus, it is difficult to have a threshold 
to distinguish between the one-class and two-class block by 
using this measure. 

B. Performance Analysis 

In this subsection, we analyze the performance of the 
QMBTC operation for the case when the input color data set is 
composed of two equiprobable 3-D normal distributions. For 
each distribution, the components of data have equal variances. 
The covariance matrix C of these distributions is given by 

(a 2 0 0 \ 
C= 0 a 2 0 
\ 0 0 a 2 ) 

where a represents variance. The probability density function 
of these normal distributions are then denoted as 

//(/) = (27r^ 2 \C\ l ^r l oxp(-£(7 - m ( yc- l {l - mi)) 
for I = 0. 1 

where J = (A./*, /a) is the vector representing the color 
value of one data point mi = (/it.Mf*Mf)» f° r i = 0, 1, are the 
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mean vectors for each distribution, and \C\ is the determinant 
of covariance matrix C According to the moment theorem for 
normal distribution [14], quaternion moments m t , 7h>. and m 3 
can be calculated as 

rh { = ± - (/*o + pi) • * + £ • (/i 0 + Mt) • j 

+ | ■(/*© + Mi)- * 
m 2 =3 (* 2 + §•(*> + /'?)) 



ff 2 -((/*o + Mi) + §- (Po + /*?))•« 



+ 5 * * " ((MO + Ml) + § • (Mo + M?)) 'J 
+ § ; <r 2 - ((mo + Ml) + § • (Mo + Mi)) • *- 

Then. i 0 and £i can be obtained by using (16) and (17): 

for ti = 1.2.3 



-o#? = -g- • (3 • cia — y0c| 3 — l2coo) 
*u = • (3 • c w + v/9c? 3 - 12coo) 



(19a) 



where 



coo =3(mo - Mi - <* 2 ) 

-(Mo + Mi)(* 2 + | - (Mo - Ml) 2 ) 
Cl3 (3 (mo -Mi) 2 ) 



(1 9b) 



Suppose that two normal distributions are well separated 
with the assumption that (p^ — Mi) ^ ^ and (mo * Mi) > <r 2 . 
From (19b), C13 is approximated as 

C13 ~ -(/i 0 +Ml). 

Besides, we understand from (19a) that the values of z lht and 
zi„ of (19a) are positive. This situation means that the input 
data set could be separated into two classes by the decision 
boundary /', which approximately passes through the point 
(((Mo + Mi)/2), ((mo + Mi)/2), ((mo + Mi)/2))). On the other 
hand, when two normal distributions are not well separated, 
there might be only one class generated by the QMBTC. 
Considering the special case when two normal distributions 
are identical, that is pty = p\ % zq and z\ become 



*o» = i • (Mo - v ^4Mg + 9or2 ) 

*u = \ • (mo + v/^fTo^) 



for v =s 1,2.3 



from which we notice that 2 0 ». is negative and the decision 
boundary /' passes through the point ((mo/3), (/i<)/3), (mo/3)). 
In this case, we can expect that QMBTC produces one 
output class with the quantized level z L being approximately 
(Mo*Mo-Mo) when a is small enough. 

In agreement with the examples in the last subsection, the 
above analysis reveals that QMBTC can decide the number 
of output classes according to the distribution of color values 
inside the pixel block. By this characteristic, QMBTC can 
achieve higher compression ratios since pixel blocks with 
similar color values need only one reproduction color, and 
the computed bit map can be replaced by a 1 bit reference 
indicating that one color clustering happened. However, the 
thresholding methods utilized by color BTC algorithms [7]-[9] 
always generate two output classes and no bit map saving. 

IV. Application to Color Image Compression 
In this section, we will introduce an approach of weighted 
color space to implement the QMBTC. To further compress 
the color image, a post- BTC data compression scheme, which 
combines thevblock merging and color codebook generation 
steps, is presented. 

A. The Weighted Color Space Approach 

As is known, the RGB space has extensive correlation 
among color components. Besides, pixels within the block are 
likely to have spatial correlation, which results in similar color 
values, except for the edge blocks. These two factors cause 
too many one-class pixel blocks by applying the QMBTC on 
the RGB space. Even though the compression ratio is high 
in this space, the picture quality of the reconstructed image 
is not acceptable from our empirical results. To alleviate this 
situation, we first transform the RGB space to the Yrg space. 



Y = 



_ (R + G + B) 

3 
R 

r (R + G + B) 
G 



9 = 



(R + G + B)' 



In this space, the color components are decorrelated into 
two parts, the achromatic and chromatic components. The 
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(0. 147, 79. 90) (0. 134. 74. 93) 
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(0.142.89,92) (0,13,14.34) 

(0,17,14,19) (0.145.93.93) 



(0,144,79,101) (0,150.86,99) \ 

(0,16,18.39) (0.147.91.106) | 

(0,15.17,38) (0.155.97,117) I 

(0, 143,92, 117) (0. 142.90. 122)/ 
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(a) 



(c) 




(b) 



id) 



Hg- I- Original testing color image., referred ,o as (a , "Lena", fb) "Peppers". 

achromatic Y component represents the intensity of the image 
: d S2 a,M P rinc jPf "- of the KL transfer^ 
the RGB space as reported by [15J. The r and g chromatic 

"s^et rcPreSem n0rmaH2ed Chr0ma,iC infonn ^ 

ulh^^IrT r" f0r , ,he associa,ed co,or com P~ 

S£ v i ~ • L ^ determ "»ed by the smoothness of 
*e Y values .ns.de the pixel block. If the sample mean value 
of Y ms.de the p.xel block is denoted as F. tj are defined a! 



t^i = 1 — exp | 



(c) "Scene", (d) "J e r. 



TABLE I 

n^ R ^ H ° LDED BLOCK DlSTR,BUTION OF 

QMBTC under Different Color Spaces 



Color 
Spaces 



One-CUu 
Blocks 



Two-Class 
Blocks 



_«p(-|) 

2 
2 

:r = Jv r E IO»-F| 



U?2 = 



LENA 



PEPPERS 



SCENCE 



JET 



I 
U 



I 

0 



I 



(20) 



11858 
9408 



12199 
9408 



10396 
7576 



11430 
9670 



4526 
6976 



4(85 

6976 



5988 
8808 



4954 
6714 



where p is a sca| mg factor md Y(n) h ^ 

nth p,xel ,n the block. The weighted color components are 

ee S tT ai T • J h,Ch ^ aPP ' ied ,0 thC QMBTC F ™ <20). we 
mocK. When T = 0, there is a constant Y value inside the 
p.xel block. The CCC algorithm f7J, which processes the }^ 
component only, would not produce satisfactory results in this 
suuat.ons.nce there might be different colors existing Within 
the pixel block. Nevertheless, the weighted color component 
approach of QMBTC can solve this problem. STS 
we understand that the weighted r and, componeTwolla 
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- P n F g when y th US,n8 QMBTC °" ima *« ^ 
» tig . 1 when the input color spaces are the RGB snac* 

and we.gh.ed Yr 3 space with „ = 3.0. respectively ^cl 

seen from Table I that the arrangemen'of «, J by 

assists .„ .mproving the situation of too many one-class pte 
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Fig. 2. Block-merging siage of ihe posi-BTC data compression scheme.. 

blocks for each test image. Besides, from empirical results, 
we noticed that the total computing time of QMBTC was 
enlarged about 6% by the weighted Yrg space approach on 
average. 

B. Post-BTC Data Compression 

To do color image compression, the proposed QMBTC is 
employed to each pixel block. After the number of classes 
has been decided, the centroid of each class is assigned as a 
reproduction color of the pixel block. The bit maps and the 
reproduction colors of all pixel blocks will then be sent to 
the receiving end and stored for decoding. However, since the 
total number of reproduction colors and bit maps for the whole 
image is still large, a post-BTC data compression scheme, 
which consists of two stages: 1) block-merging stage, and 2) 
color codebook generation stage, is presented in this paper to 
further compress the image. 

In the block-merging stage as illustrated in Fig. 2, the input 
image will be first divided into nonoverlapping groups of 16 
pixel blocks. In each group of pixel blocks, all blocks will 
be examined to see whether they can be merged as a larger 
pixel blocks, called a macroblock. If every block in the group 
is a one-class block, this group of blocks will be marked 
as a one-class macroblock. Otherwise, this group of blocks 
cannot be merged. For the case of 4x4 pixel block, the size 
of the merged macroblock becomes 16 x 16. For a merged 
macroblock, we then extract the average color value of each 
pixel block to form a 4x4 size block. On this block, QMBTC 
is executed to produce one or two reproduction colors for the 
entire merged macroblock. When one reproduction color is 
generated, the whole merged macroblock can be referenced 
with only one bit. Otherwise, when two reproduction colors 
are generated, each pixel block within a merged macroblock 
chooses its reproduction color according to its average color 
value belonging to the resultant color clustering. 

Concerning the color codebook generation stage, a color 
quantization method proposed by Pei and Cheng [16], which 
partitions the color space of the image sequentially according 
to the order of color components, is first employed. Then, 
for each reproduction color generated by the QMBTC, its 
associated color index is generated by finding the nearest 
neighbor entry of the color codebook. The proposed post-BTC 
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Fig. 3. Process of the post-BTC data compression scheme. 

data compression process is shown in Fig. 3. To decode each 
compressed pixel block, the bit map, color indexes of each 
pixel block, and the color codebook are the only required 
elements. The receiver can display the decoded data through 
a lookup table manner in real time. 

V. Experimental Results 
In this section, we conducted the experiments on im- 
ages of Fig. 1, called "Lena," "Peppers," "Scene," and "Jet." 
These test color images with 512 x 512 size are coded 
at 8 bits/component. The average peak signal-to-noise-ratio 
(APSNR) criterion described as follows is used to judge the 
picture quality. 



APSNR = 10 log 



3 » N * 255 2 
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Fig. 4. (a) Output bit-rate distribution of the QMBTC versus the scaling factor p for four test 
versus the scaling factor p for four test images. 



images, (b) Output APSNR distribution of the QMBTC 



III shows, the bit rate can be further decreased to 1 13 ners" "Seen*" nnri «w *~ u ^ A , L 

bits/pixe, on average. In addition, from empirical result* 3 Sme'blocTng exist on imaged 

Z^Z^T^V 0 ^ T£T ,imC ° f QMBTC WaS " Lena " and " Pte PP«.- ** qualms accede to 

enlarged about 3% by post-BTC data compression on aver- the human eye at a distance of four-five times the picmre 

age. The associated reconstructed images for "Lena". "Pep- height picture 



DOCICt <XP BS9448A I > 



> 



II-KK TRANSACTIONS ON COMMUNICATIONS, VOL. 45. NO. V MAY |W7, 




(a) (b) 




<ei 

Fig. 5. Original and reconstructed images of "Lena." (a) Original image, lb) Reconstructed image by the QMBTC with block size 4x4. (c) Thresholded 
block distribution by the QMBTC. Id) Reconstructed image by the CCC algorithm with block size 4x4. le) Reconstructed image by the Kurita and 
Otsu algorithm with block size 4x4. 
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fig. 6. Original and reconstructed images of "Jet." (at Original image. <b) Reconstructed image by the QMBTC with block size 4x4. (c) Thresholded 
block distribution by the QMBTC. (d) Reconstructed image by the CCC algorithm with block size 4x4. (e) Reconstructed image by the Kurita and 
Olsu algorithm with block size 4x4. 
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Fig. 7. Reconstructed images of QMBTC after post-BTC data compression, (a) Reconstructed image of "Lena." (b) Reconstructed imaee of "Peppers 
(ci Reconstructed image of ■'Scene." (d) reconstructed image of "Jet." rr 



VI. Conclusions 

In this paper, we present a color BTC algorithm, called the 
QMBTC, to perform color image compression. The proposed 
algorithm can produce one-class or two-class pixel blocks 
through quaternion-moment-preserving thresholding. From the 
derived formulas, we know that the computation complexity of 
the QMBTC is of order iV, the size of pixel block. Empirical 
results, which compare the APSNR's and bit rates of the 
proposed algorithm with those of existing single bit map color 
BTC algorithms, show that a 30% saving is achieved on 
average by the proposed algorithm, while no significant degra- 
dation of picture quality is perceived. To further compress 
color image, a post-BTC data compression scheme, which 
combines the block-merging and color codebook generation 
steps, is employed. The bit rate is then dropped to 1.13 
bits/pixel on average with acceptable output picture quality. 
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